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The interplay between magnetism and metal-insulator transitions is fundamental to the rich 
physics of the single band fermion Hubbard model. Recent progress in experiments on trapped 
ultra-cold atoms have made possible the exploration of similar effects in the boson Hubbard model 
(BHM). We report on Quantum Monte Carlo (QMC) simulations of the spin-1 BHM in the ground 
state. For antiferromagnetic interactions, {U2 > 0), which favor singlet formation within the Mott 
insulator lobes, we present exact numerical evidence that the superfluid-insulator phase transition 
is first (second) order depending on whether the Mott lobe is even (odd). Inside even Mott lobes, 
we examine the possibility of nematic-to-singlet first order transitions. In the ferromagnetic case 
{U2 < 0), the transitions are all continuous. We map the phase diagram for U2 < and demonstrate 
the existence of the ferromagnetic superfluid. We also compare the QMC phase diagram with a third 
order perturbation calculation. 

PACS numbers: 05.30.Jp, 03.75.Hh, 67.40.Kh, 75.10.Jm 03.75.Mn 



The single band fermion Hubbard model (FHM) offers 
one of the most fundamental descriptions of the physics 
of strongly correlated electrons in the solid state. The 
spinful nature of the fermions is central to the wide 
range of phenomena it displays such as interplay be- 
tween its magnetic and transport propcrtiesji Such com- 
plex interplay is absent in the superfluid to Mott insula- 
tor transition^ii^i^i^ in the spin-0 Boson Hubbard model 
(BHM). However, purely optical traps^ can now confine 
alkali atoms ^'^Na, '^^K, and ^^Rb, which have hyper- 
fine spin F = I, without freezing Fz- As in the fermion 
case, the nature of the supcrfiuid-Mott insulator (SF-MI) 
transition is modified by the spin fluctuations which arc 
now allowed. Initial theoretical work employed contin- 
uum, effective low-energy Hamiltonians and determined 
the magnetic properties and excitations of the superfluid 
phases ii 

To capture the SF-MI transition it is necessary to con- 
sider the spin-1 Bosonic Hubbard Hamiltonian. 

H= -t "^{al^aja + h.c) + ^"^h,{h,~l) 
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The boson creation (destruction) operators {dia) have 
site i and spin a indices, a = 1, 0, —1. The first term de- 
scribes nearest-neighbor, (ij), jumps. The hybridization 
t = I sets the energy scale, and wc study the one di- 
mensional case. The number operator fii = hi^ ~ 
So- ^la^iiy counts the total boson density on site i. The 
on-site repulsion C/q favors states with uniform occupa- 
tion and competition between C/q and t drives the MI-SF 
transition. The spin operator Fi = „i cil^F^^iaia' , 

with F^ra' the standard spin-1 matrices, contains fur- 
ther density-density interactions and also interconver- 
sion terms between the spin species. We treat the sys- 



tem in the canonical ensemble where the total particle 
number is fixed and the chemical potential is calculated, 
^i{N) = E{N+1)-E{N) where eIn) is the ground state 
energy with N particles. This holds only when each en- 
ergy used is that of a single thermodynamic phase and 
not that of a mixture of coexisting phases as happens 
with first order transitions. It is therefore incorrect to de- 
termine first order phase boundaries with the naive use of 
this method. In the present case there is the added sub- 
tlety that the three species are interconvertible. These 
issues will be addressed in more detail below. 

Several important aspects of the spin-1 BHM are re- 
vealed by analysing the independent site limit, t/Uo ~ 0. 
The Mott-1 state with rii ~ \ on each site has £m(1) = 0. 
In the Mott-2 state with rii = 2, the energy is fM(2) = 
Uq — 2U2, if the bosons form a singlet, F = 0, and is 
£m2 = Uo + U2 if F ^ 2. Thus U2 > favors singlet 
phases while U2 < favors (on-site) ferromagnetism. 
This applies to all higher lobes as well. In the canon- 
ical ensemble, the chemical potential at which the sys- 
tem goes from the nth to the {n + l)th Mott lobe is 
/i(n n + 1) = £n+i - £n- 

First consider U2 > The energy of Mott lobes at 
odd filling, rio, is f A/(no) = UoUoirio - l)/2 -I- 172(1 - Uo) 
while at even filling, jic, £j\/(no) = Uonc(nc — l)/2 — ncU2- 
Therefore, the boundaries of the lobes, going from lower 
to higher filling, are /i(rio (ne+l)) = nJjQ and ^{rio 
(rio + 1)) = noUQ — 21/2- This demarks the positions of 
the 'bases' of the Mott lobes in the {t/Ua, fi/Uo) ground 
state phase diagram. For U2 > the even Mott lobes 
grow at the expense of the odd ones, which disappear 
entirely for Uq = 2U2- 

For U2 < 0, the ground state is ferromagnetic (max- 
imal F: F^ = n{n + 1)) which gives for all Mott lobes 
Suin) = n{n — l)(t/o + C^2)/2. Consequently, the bound- 
ary of the nth and (n + l)th Mott lobes /i(n — > n -I- 1) = 
n([/o + U2)- The bases of both the odd and even Mott 
lobes shrink with increasing IC/2I, in contrast to the 



U2 > case where the even Mott lobes expand. 

Mean-field treatments of the lattice model capture the 
SF-MI transition as the hopping t is turned on, and have 
been performed both at zero and finite temperature 
Even when U2 ~ 0, the spin degeneracy alters the nature 
of the transition. For U2 > 0, the order of the phase 
transition depends on whether the Mott lobe is even or 
odd. These mean field calculations assume a non-zero 
order parameter (aia-), which cannot be appropriate in 
d=l or in d=2 at finite T. Therefore it is important 
to verify these predictions for the qualitative aspects of 
the phase diagram, especially in low dimension. A quan- 
titative determination of the phase boundaries requires 
numerical treatments. Indeed, DMRG^i and Quantum 
Monte Carlo (QMC)i^ results for C/2 > and d=l re- 
ported the critical coupling strength and showed the odd 
Mott lobes are characterized by a dimerized phase which 
breaks translation symmetry. 

For U2 < 0, the nature of the SF-MI transition does 
not depend on the order of the Mott lobe while for U2 > 
it is predicted to be continuous (discontinuous) into 
odd (even) lobes. Consequently, it suffices to study the 
first two Mott lobes both for U2 > and J72 < to 
demonstrate the behavior for all lobes. Furthermore, in 
what follows we will focus on the case \U2/Uo\ = 0.1 in 
order to compare our results for U2 > with Rizzi et aJii 

Here, we will use an exact QMC approach, the SGF 
algorithm with directed update, to study the spin-1 BHM 
in d=l for both positive and negative 1/21^ 

For U2 > (e.g. ^^Na), which favors low total spin 
states. Fig. [T] shows the total number density p = N/L 
against the chemical potential, fx, for Uq = lOt and 
U2 = t. It displays clearly the first two incompressible 
MI phases. In agreement with the t/Uo = analysis, U2 
causes an expansion of the second Mott lobe, p = 2, at 
the expense of the first, p = 1. Our Mott gaps agree 
with DMRG result&ii to within symbol size. However, 
the p versus /i curve in Fig. [T] does not betray any ev- 
idence of the different natures of the phase transitions 
into the first and second Mott lobes. In particular, for a 
spin-0 BHM, first order transitions are clearly exposed by 
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the appearence of negative compressibility 
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dp/dp, 



which is not present here. The transition into the sec- 
ond Mott lobe is expected to be first order and driven by 
the formation of bound pairs of bosons in singlet states. 
Therefore, in the canonical ensemble, we expect that near 
this transition there will be phase coexistence between 
singlets arranged in a Mott region and superfluid. 

The nature of the transitions is revealed by the evo- 
lution of the spin populations in the system as p in- 
creases. Since the singlet wavefunction is |0, 0) = 

y273|l,l)|l,-l) - yT73|l'0>|l'0>. tliis state has p+ = 
po = P- • We plot in the inset of Fig. [1] the population 
fractions, Nq/N and N^/N = N+/N versus the total 
density. We see that as p increases, N+/N and Nq/N 
oscillate: When N is even, singlet bound states of two 
particles try to form drawing the values of N^/N and 
Nq/N closer together. However, singlets form fully, mak- 
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FIG. 1: (Color online) p{fi) exhibits Mott plateaux: gapped, 
insulating phases at commensurate fillings. Inset: N+/N and 
iVo /N vs p = N/L. Singlets form where N+ /N and A^'o /N are 
equal (see text). 



ing N+ = Nq, only close to the second Mott lobe, p = 2, 
where we clearly sec N^/N = Nq/N for a range of even 
values of N. On the other hand, when N is odd, singlets 
cannot form and the spin populations are much farther 
from that given by the singlet wavefunction. In the ther- 
modynamic limit and fixed N, one expects true phase 
separation into p = 2 singlet MI regions and p < 2 SF re- 
gions. For a finite system, phase separation commences 
for (even) fillings where we first have N+/N = Nq/N, 
i.e. 1.5 < p < 2 and similar behavior for p > 2. An- 
other interesting feature in the inset of Fig. [1] is that the 
difference between iV^/A^ and Nq/N, for even N, de- 
creases linearly as the density approaches the transition 
at p = 1.5. No such behavior is seen as the first Mott 
lobe is entered from below or above p = I', the transition 
is coninuous as predicted for odd lobes. 

The boxes in Fig. [T]show the values of p corresponding 
to phase coexistence rather than to one stable thermo- 
dynamic phase. It is, therefore, clear that the canonical 
calculation of the phase boundaries, i.e. simply adding a 
particle to, or removing it from, the MI, is not applicable 
in the presence of a first order transition. 

Figure [1] reveals the SF-MI transition at fixed, suffi- 
ciently large Uq when p is varied. In Fig. [2] we show the 
transition when the density is fixed, p = 2, and Uq is 
the control parameter with U2/UQ fixed at 0.1 (main fig- 
ure) and 0.01 (inset). Singlet formation is clearly shown 
by (F^) 0, as t/UQ decreases and the second MI is 
enteredfi^ Indeed, the origin of the first order transition 
into even Mott lobes, as the filling is tuned, is linked 
to the additional stabilization of the Mott lobe associ- 
ated with this singlet energyJ^ The superfluid density, 
Ps = L{W'^) /2tp, where W is the winding number, is a 
topological quantity and truly characterizes the SF-MI 
phase transition which is continuous in this case. As L 
is increased from 10 to 16 and 20, the vanishing of ps 
gets sharper. We find that the critical value of t/Uo for 
the p = 2 lobe is somewhat less than that reported in 
DMRGii indicated by the dashed line. We believe this 
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FIG. 2: (Color online) The square of the local moment, {F'^), 
and the superfluid density, ps, vs t/Uo for p = 2. In the 
Mott lobe, (F^) signaling singlet formation and ps 
indicating an insulator. The dashed line indicates the critical 
t/Uo from Rizzi et aS^. Inset: Same but with U2/U0 = 0.01. 



is because with DMRG the phase boundaries were ob- 
tained using finite differences of the energy with small 
doping above and below commensurate filling. As dis- 
cussed above, this is not appropriate for a first order 
transition. 

For 2dU2/Uo < 0.1, and d = 2,3, mean fieldi^^i^ pre- 
dicts that, when t/U§^ 



^/U^JIdU^ is in the MI, then 
Mott lobes of even order are comprised of two phases: 
(a) the singlet phase for t/Uo ^ V^o^ ^^'^ (b) a nematic 
phase for t/U§'^ < t/Uo < t/U;^, where t/U§ is the tip 
of the Mott lobe. Inside the lobe, the nematic-to-singlet 
transition is predicted to be first order which raises the 
question: are the singlet-to-SF and the nematic-to-SF 
transitions of the same order? Figure [2] shows that the 
SF-MI transition, ps 0, occurs at larger t/Uo than 
singlet formation, {F'^) — > 0, both for U2/U0 = 0.1 and 
0.01. The passage of {F'^) to zero gets sharper for smaller 
U2/U0 but remains continuous, not exhibiting any signs 
of a first order transition. We have verified this for 
U2/U0 = 0.1, 0.05, 0.01, 0.005. Furthermore, the insen- 
sitivity of {F'^) to finite size effects indicates that it does 
not undergo a continuous phase transition. We have also 
verified that for the second Mott lobe, the SF-MI transi- 
tion is first order regardless whether t/Uo is less or greater 
than t/U§^. We conclude that while ps ^ is a contin- 
uous critical transition, {F'^) — > is a crossover not a 
phase transition. This, of course, does not preclude the 
possibility of a first order transition for d = 2,3. 

Whereas ^-^Na has positive U2: ^''Rb has U2 < 0, 
leading to different behavior. We begin with p ver- 
sus p in Fig. [3] Unlike the U2 > case, the SF- 
MI transitions are continuous for both even and odd 
Mott lobes: the inset shows that the spin populations 
do not oscillate as for U2 > 0- The population ra- 
tio, po = 2p+ can be understood as follows. As shown 
above for t/Uo ^ 0, maximum spin states are favored 
when U2 < 0. So, when a site is doubly occupied, the 
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FIG. 3: (Color online) The density, p versus chemical poten- 
tial, p, exhibits the usual Mott plateaux at commensurate 
filling for U2 < 0. Inset: The spin population fractions show- 
ing Po ~ 2p+. They do not exhibit oscillations like in Fig. [1] 



spin-2 state is favored. But, since our study is in the 
^totai „ _ — sector, the wavefunction of the 
spin-2 state is |2,0) = l/\/3|l, 1) |1, -1)-F 0) |1, 0) 

and thus po = 2p+ = 2p_. 

As discussed above, U2 < favors 'local ferromag- 
netism', namely high spin states on each of the individ- 
ual lattice sites. As with the FHM, the kinetic energy 
gives rise to second order splitting which lifts the de- 
generacy between commensurate filling strong coupling 
states with different intersite spin arrangements. We can 
therefore ask whether the local moments order from site 
to site: Do the Mott and superfluid phases exhibit global 
ferromagnetismi^? To this end, we measure the magnetic 
structure factor. 



(2) 



where a — x or z. Figure |4] shows S-xxiq) in the super- 
fluid phase at half-fiUingji^ The peak at g = grows 
linearly with lattice size, indicating the superfluid phase 
does indeed possess long range ferromagnetic order. We 
find that the MI phase is also ferromagnetic. 

To determine the phase diagram, we scan the density 
as in Fig. [3] for many values of Uo with U2/U0 constant 
(—0.1 in our case). The resulting phase diagram is shown 
in Fig. [5l Comparison of data for two lattice sizes demon- 
strates that finite size effects are small. 

Early in the evaluation of the phase boundaries of 
the spin-0 BHM it was observed that a perturbation 
calculation^ agreed remarkably well with QMC results^ 
We now generalize the spin-0 perturbation theory to spin- 
1 and show a similar level of agreement with the QMC 
results. If we assume the system always to be perfectly 
magnetized, then n bosons on a site will yield the largest 
possible spin, F^ = n{n + 1). Consequently, the in- 
teraction term in the Hamiltonian, Eq.((T]), reduces to 
(Uq — U2) ~ givmg a Hamiltonian identi- 

cal to the spin-0 BHM but with the interaction shifted 
to {Uo — U2)/2. One can then repeat the perturbation 
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FIG. 4; (Color online) The magnetic structure factor 5'xx(g) 
for U2 < and p = 0.5 exhibits a sharp g = peak indicating 
ferromagnetic order. 
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FIG. 5: (Color online) Phase diagram of the spin-1 BHM for 
U2/U0 = -0.1. 



expansion to third order in t/ {Uq — U2) to determine the 



phase diagram The result is shown as the dashed Une in 
Fig. [Hand is seen to be in excellent agreement with QMC. 
The agreement further suggests that the finite lattice ef- 
fects in the phase diagram are small. Such a perturbation 
calculation is not possible for the [/q > case since 
depends on the phase, SF vs MI, and on the order of the 
MI lobe. 

The (dipolar) interactions between spinful bosonic 
atoms confined to a single trap have been shown to give 
rise to fascinating "spin textures" .— An additional opti- 
cal lattice causes a further enhancement of interactions, 
and opens the prospect for the observation of the rich 
behavior associated with Mott and magnetic transitions, 
and comparisons with analogous properties of strongly 
correlated solids . ^^'^^ Here, we have quantified these phe- 
nomena in the one-dimensional spin-1 BHM with exact 
QMC methods. We have shown that, for U2 > 0, the MI 
phase is characterized by singlet formation clearly seen 
for even Mott lobes where (F^) — > as fJo increases. 
We also showed that the transition into odd lobes is con- 
tinuous while that into even lobes is discontinuous (first 
order). We emphasized that the naive canonical determi- 
nation of the phase boundaries is not appropriate for a 
first order transition. For U2 < 0, we showed that all MI- 
SF transitions are continuous and that both the SF and 
MI phases are ferromagnetic. The phase diagram in the 
{fi/Uo,t/Uo) plane obtained by QMC can be described 
very accurately using third order perturbation theory. 
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